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^ : ABSTRACT 

| We study issues concerning perturbative integrability of 5\£ = 6 Chern-Simons theory at planar 
and weak 't Hooft coupling regime. By Feynman diagrammatic s, we derive so called maximal- 
ranged interactions in the quantum dilatation generator, originating from homogeneous and 
inhomogeneous diagrams. These diagrams require proper regularization of not only ultraviolet 



but also infrared divergences. We first consider standard operator mixing method. We show that 
homogeneous diagrams are obtainable by recursive method to all orders. The method, however, 
is not easily extendable to inhomogeneous diagrams. We thus consider two-point function 
method and study both operator contents and spectrum of the quantum dilatation generator up to 
six loop orders. We show that, of two possible classes of operators, only one linear combination 
actually contributes. Curiously, this is exactly the same combination as in =4 super Yang- 
Mills theory. We then study spectrum of anomalous dimension up to six loops. We find that 
the spectrum agrees perfectly with the prediction based on quantum integrability. In evaluating 
the six loop diagrams, we utilized remarkable integer-relation algorithm (PSLQ) developed by 
Ferguson, Baily and Arno. 



1 Introduction 



The AdS/CFT correspondenceQ]| continues revealing remarkable relations between gauge and 
gravity theories. The most extensively studied so far is the correspondence between the four- 
dimensional 3\£ = 4 superconformal Yang-Mills theory and the Type IIB superstring theory on 
AdSs x § 5 [1J. Importantly, both theories admit Lagrangian formulation, which involve two 
coupling parameters: the rank N of the gauge group G and the 't Hooft coupling constant X 
for the former, and the string coupling g s and the curvature scale R (as measured in string unit) 
for the latter. As such, perturbatively, one can compute physical observables in both theories 
in double series of the respective parameters and test the correspondence by comparing a given 
observable extracted from each sides. In the planar limit, N — > °° and g s — > 0, remarkable 
agreement was discovered between the two sides for a variety of observables. The agreement is 
largely attributed to the integrability structure B- lfTTl . 

Recently, a Type IIA counterpart was discovered and added to the list of the AdS/CFT 
correspondence. The ABJ(M) theory is (2+l)-dimensional fA£ = 6 superconformal Chern- 
Simons theory and was proposed as the holographic dual to the Type IIA superstring theory 
on AdS4 x CP 3 lfT8l[T9ll . A question of interest is whether the two sides in this Type IIA coun- 
terpart also have an integrability structure. Recently, positive indications toward the quantum 
integrability were accumulated ll2"0l - 113711 . At strong coupling side, Lax pair construction of the 
integrability was shown at leading order[20J. At weak coupling side, there were more indica- 
tions. At two loops, spin chain Hamiltonian was computed explicitly for the SO(6) subsector 
and was shown integrable lETl 1221 1231 . See also Ref. [27]. At four loops, spectrum of the spin 
chain Hamiltonian was shown to agree with the prediction of the integrability [|35l l36ll . 

In this paper, we further continue our previous investigations [|22l l25l l35l concerning inte- 
grability of the 9l = 6 Chern-Simons theory in the weak 't Hooft coupling regime by comput- 
ing the dilatation operator to six loops. Given that the integrability is in place up to four loops, 
why bother six loops? We contend that there are two important aspects that arise beginning at 
six loops and beyond: operator contents and recursive structure of the long-range spin chain. 
With these two issues on focus, we shall test the integrability of the 3\£ = 6 Chern-Simons 
theory at six loop order. As in [|35ll . we shall focus on magnon excitation in the SU(2) sub- 
sector, compute operator structure and spectrum of the spin chain Hamiltonian up to six loops 
and check them against the prediction based on the integrability and the centrally extended 
[psu(2|2) ©psu(2|2)] x M 2 ' 1 superalgebra of excitation symmetry. 

The off-shell psu(2|2) superalgebra of the excitation symmetry is spanned by the two su(2) 
rotation generators £ a p, the supersymmetry generator 0™ and the superconformal gener- 
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ator ©a- The off-shell configuration is characterized by s/(2,]R) central charges .ft* [fT5l . 
Their (anti)commutators are given by [fT5l 

m, r] =v b r- \m c , [£V 3 Y ] = sjr - 

{0« nj} = e a h ah &, {6S, ©{} =e a p£ a ^. (1.1) 

The central charges £ is related to the energy by E = £, while .ft, .ft* introduced at off-shell are 
related to momentum of the magnon. Acting on a magnon transforming in the fundamental 
representation, closure of the superalgebra leads to the relation among the central charges 

E 2 = <£ 2 = ~+4SUi*. (1.2) 

The central charges .ft, .ft* are in turn related to an excitation momentum P. More generally, a 
bound-state of Q elementary magnons transforming in higher-dimensional representations can 
be studied. The off-shell analysis was sufficient to determine the dispersion relation. It is 

E(P) = l -^Q 2 + \6h\\) sin 2 1, (1.3) 

where h(X) is a function of the 't Hooft coupling parameter X. 

Restricting to large N limit and 50(6) sector of the OSp(6\4,W) superconformal symme- 
try group, the quantum dilatation operator was computed explicitly at two loops from which 
an integrable alternating spin chain Hamiltonian and Bethe ansatz equations were identified 
ET1I221I251 . Aspects of the integrability were explored further beyond two loops by focusing 
on diagrams generating maximal-ranged terms. These are the diagrams in which interaction 
vertices range over lattice sites of the spin chain Hamiltonian maximally. In [|35ll . we computed 
four-loop contribution to these terms and found that the spectrum fits with the prediction based 
on the integrability and the excitation symmetry. 

At each order in perturbation theory, depending on the range the 'spin' flavors at different 
sites are permuted, maximal-ranged terms in the dilatation operator are further classifiable into 
maximal-shuffling and next-to-maximal-shuffling terms. In 91 — 4 SYM theory, it was found 
by Gross, Mikhailov and Roiban [3] that maximal- shuffling terms are computable recursively. 
Inspection of relevant Feynman diagrams indicates that all diagrams contributing to maximal- 
shuffling terms are related by a recursion relation and hence resummable to an exact all-order 
result. On the other hand, diagrams contributing to nonmaximal-shuffiing terms are combinato- 
rially too complicated to be resummable. One might anticipate that a similar argument exists for 



2 



the $C =6 ABJ(M) theory since conformal interactions are tightly constrained by large amount 
of supersymmetry. Indeed, we shall find that the maximal- shuffling terms in the dilatation oper- 
ator originates from the same class of skeleton diagrams which we call homogeneous diagrams. 
We were able to perform all-order resummation and show that they match with the structure 
of dilatation operator predicted by integrability. The nonmaximal-shuffiing terms receive con- 
tribution from another class of skeleton diagrams which we call inhomogeneous diagrams. As 
these diagrams are not recursively resummable and afflicted with potential infrared divergences, 
we need to resort to an alternative approach for direct evaluation. In the second half of work, 
we thus adopt two-point function method and compute six-loop contribution to maximal-ranged 
interactions from both homogeneous and inhomogeneous diagrams. This method amounts in 
dual Type IIA string theory to deriving time evolution Hamiltonian of a single non-interacting 
string. 

Key results of our work point to the following s. The dilatation operator at six loops are 
extractable free of infrared divergences from two-point correlation functions of gauge invariant 
operators, which was already utilized in our earlier study at four loops 0511 . Moreover, the 
maximal-ranged interactions are consistent with the integrability and that, curiously, operator 
contents of the long-range spin chain Hamiltonian is identical to those of the 3\£ = 4 super 
Yang-Mills theory, viz. Inozemtsev spin chain Il38l l8ll. 

This paper is organized as follows. We begin with description of the expected spectrum 
based on the integrability and prediction for the maximal shuffling coefficients to all orders. 
In section 3, we recapitulate all-loop computation of maximal shuffling terms in 5\£ = 4 super 
Yang-Mills theory, obtained first by Gross, Mikhailov and Roiban in J3). We then extend the 
method to 5\£ = 6 Chern-Simons theory and find two classes of diagrams contributing to the 
maximal shuffling. The first class of diagrams, homogeneous diagrams, is computable by a 
straightforward extension of the method in [0 and yields a result exactly parallel to the 5\£ = 4 
super Yang-Mills theory. The second class, inhomogeneous diagrams, is not computable by 
the method of or variants of it: these diagrams are afflicted with infrared divergences. We 
conclude that all-order derivation for the maximal shuffling part of the dilatation operator is not 
straightforward in 5\£ = 6 Chern-Simons theory. We thus resort to computing operator contents 
and spectrum of their anomalous dimensions order by order in perturbation theory. In section 4, 
we study 6-loop contribution to the anomalous dimension directly defined by two point correla- 
tion functions of operators. In section 5, we identify the matrix structures of the homogeneous 
and the inhomogeneous maximal-ranged interactions at 6-loops. In section 6, we compute the 
maximal-ranged part of the 6-loop Hamiltonian and confirm the prediction based on the inte- 
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grability. We also show that only one particular choice arises among two possible maximal 
shuffling operators. This operator is the same as the one arising in the 5\£ = 4 super Yang-Mills 
theory. In section 7, we extend our results to parity non-conserving 5\£ = 6 ABJ theory, whose 
gauge group is U(M) x U(iV) with N < M < N + k. The last section is devoted to the conclud- 
ing remarks. In the appendices, we relegate several technical details. Appendix A illustrate 
a comparative calculation of inhomogeneous term in operator-mixing method. Appendix B 
presents several lattice operator identities. Appendix C discusses derivation of skeleton 4-loop 
diagrams. Appendix D explains implementation of numerical integration and the remarkable 
integer relation PSLQ algorithm. 

2 Hamiltonian and Spectrum From Integrability 

We begin with consequence of the quantum integrability and the off- shell superalgebra to the 
spin chain Hamiltonian of a sub-sector of our interest. Consider single-trace bosonic operators 
in the 91 = 6 ABJM theory of the type 

«P[/l/2Ws-/2i] = tr[y; 1 ^7^/ J -F 7 ! l ]. (2.1) 

Here, 2L is the total number of the sites. We shall take the infinite volume limit L — » °° and 
view a particular ordering of the operator as a lattice spin chain wave function *F. Gauge 
invariant operators place at odd sites the Yj (I, J = 1,2,3,4) elementary scalar fields trans- 
forming as 4 under the SU (4) R-symmetry and at even sites the Yj conjugate fields. Under 
the SU (4) ~ SO(6) R-symmetry, these fields transform as 4 and 4, respectively. We denote 
{Yi,Y2,Y$,Y4} as {A\,A2,B[,bI} where A a ,Ba (a,d = 1,2) transform under the SU{2)a and 
the SU(2)g subgroups of SU(4). We then consider a subset of the operators (|2.1I) . where we 
put only A a I B A fields at the odd /even sites, respectively. Explicitly, they are the following 
restricted set of operators 

¥[aia2fl304«5 "-an] =tr[A ai B &2 A a3 Ba 4 A as ■■■B dlL \. (2.2) 

Since the only possible interaction between Yj and Y' J fields is the contraction dj group theo- 
retically, there cannot be any interaction between 4 and 4 representations acting on the above 
type of states. Therefore, for this restricted set of operators, the odd- site chain and even- 
site chain are decoupled from each other. Thus, there will be A-type magnon and 5-type 
magnon propagating independently without any interactions between them. For unrestricted 
operators, there are interactions between them, but the above choice avoids unnecessary com- 
plexity in investigating the integrability. For the reference vacuum, we take the ferromagnetic 
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state |0) = *P[1 1 11 1 ••■ 1] = tr[AiBiA\BiA\ ■■■B\]. In J35J, we explained how elementary 
excitations above the reference vacuum are organized by the centrally extended \psu(2\2) © 
psu(2\2)] k M 2 ' 1 superalgebra. The first psu(2\2) acts on the flavors of the magnons formed 
by exciting odd sites of the spin chain, while the second psu{2\2) acts on the magnons at even 
sites. The bosonic su{2) subalgebra of psu(2\2) superalgebra corresponds to exciting Y2, Y4 for 
the odd sites and Y£, Yj for the even sites. In [|35l . we carried out careful study of the rep- 
resentations of the centrally extended superalgebra for the asymptotic spin chain where 2L is 
sent to infinity. The quantum integrability of the dilatation operator implies the factorization 
of multi-magnon S-matrices into product of two-magnon S-matrices satisfying the Yang-Baxter 
equations Q. For both the AB JM and the ABJ theories, we also showed ll22l 1231 that dynamics 
of the magnons on even sites and on odd sites are governed by two separate transfer matri- 
ces x alt (w,y),x alt (v, — y) for arbitrary spectral parameters M,v,yand that, using the Yang-Baxter 
equations, they are mutually commuting 



[x^, y),x alt (v, -y)] = 0, [f>, y),x alt (v, -y)] = 0, [x alt ( a ,y),x alt (v, -y)] = . (2.3) 



.alt/ 



fealt/ 



=alt 



alt/ 



^alt/ 



Their moments are 



Q„ = ar 1 lnx alt ( M ,Y) 



alt/ 



»=0 



and 



ar J inx alt (v, 



-y) 



v=0 



1,2, 



•) (2.4) 



of which Q2 + Q.2 is proportional to the dilatation operator. They all depend on the spectral 
parameter y - a nonzero value of y is an indication that the dilatation operator and all other 
higher moments are not invariant under lattice parity transformation. We found in Il22l 1231 
that y = not only for the ABJM theory but surprisingly also for the ABJ theory, which is 
parity non-conserving. From (12.41) . it followed that there are two sets of mutually commuting, 
infinitely many conserved charges 



[Qm, On] = 0, [a m , Qj = 0, [Cbn, Oj = . 



(2.5) 



It was argued 113711 that these mutually commuting conserved charges are responsible for reflec- 
tionless property of the S-matrix elements between a magnon on even sites and a magnon on 
odd sites. 

Invariance of the S-matrices under the off-shell superalgebra transformations was crucial to 
fix the representation as well as the spectrum of elementary magnons. The analysis (recapitu- 
lated in the previous section) shows that the dispersion relation of an elementary magnon takes 



'The all-loop proposal of the Bethe ansatz and S-matrices of 91 = 6 ABJM theory was put forward in Ref. Il23l . 
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the form (11.31) with Q = I and P is pseudo-momentum of the magnon and h(X) is an interpo- 
lating function of the 't Hooft coupling X. In addition, the pseudo-momentum P that specifies 
the central charge and the magnon spectrum can also be a function of the lattice momentum 
p defined by translation in the spin chain. The functional form of the interpolating function 
h(X) and the pseudo-momentum P(p) are not determinable by the symmetry alone and require 
extra inputs of explicit computations either of 3\£ = 6 Chern-Simons theory at weak coupling 
or of string worldsheet sigma-model at strong coupling. In the previous work 11351 . we found 
that P(p) = p holds up to four-loop order. In this work, we shall assume this as an input and 
proceed for computation of six-loops and beyond. 

Perturbatively, the interpolating function h 2 (k) is expandable as 

h 2 (X)=l 2 n + J^h 2e l 2i \ = X 2 (\+h 2 X 2 + h 4 X 4 + ---) (2.6) 

where, for the leading term, we use the result of the two-loop computation in ETll22~l . Recently, 
in 11361 . the next coefficient h 2 was computed to be h 2 = 4^(2) — 16. Thus, in terms of the lattice 
momentum p, the magnon dispersion relation can be expanded as 

E(P) = L^I^sin 2 ^ 

n=0 1=0 1 

= (~) + (4 sin 2 P -)x 2 + (4h 2 sin 2 | - 16 sin 4 |) A 4 

+ [e 6> 2 sin - + ^6,4 sin - + e 6)6 sm -JX +■■■ . (2.7) 

Note that eo,o = 1/2 corresponds to the classical scaling dimension of the elementary scalar 
fields Y 2 ,Y4, Y 2 \yJ and e 2n ,o = (n > 0) is required by one-third of the %t = 6 supersymmetry 
preserved by the ferromagnetic vacuum state. □ The two-loop coefficient e 2 ^ 2 = 4 was found by 
explicit computation in ETl l22l . The four-loop coefficient 64.4 = — 16 was computed in 11351 . 
while 64,2 = 4/i2 was argued in [l36l . 



2 In Chern-Simons theory, the choice of regularization method is known to be a subtle issue. Detailed study in 
utilized the dimensional reduction and obtained e2,o = at two loops. This confirms that the dimension re- 
duction is a gauge invariant and supersymmetric regulator at least up to two loop order. Whether the corresponding 
Ward identities are satisfied at higher loops is an open problem that needs to be checked. All higher loop diagrams 
involving gauge and ghost fields are afflicted by the problem but diagrams involving matter fields only are not. Our 
previous ll35l and the present works deal only with Feynman diagrams of the latter type and hence are free from 
this open problem. 
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On the other hand the expected spectrum based on the quantum integrability and the off- 
shell superalgebra representation theory is expandable as 



E(p) 



1^1 + 16/^) sin 2 1 
I + f 4 sin 2 ^U 2 +f4sin 2 ^ 



4sin 2 ^ 



+ 4 sin 



,2P 



hd — 8/z? sin — (-32 sin 

2 2 J 



2. 



(2.8) 



where we use the expansion for h 2 (k) in (12.61) . Note that the coefficients e2 n ,2n arQ fixed com- 
pletely by the dispersion relation: 



&2n.2n 



4 » J 2 "" 2 ) 1 



(2.9) 



l)\n\ 

- i P) n /A n . Therefore, at order l 2n of the 



It is the coefficient of the term (sin 2 § ) n = (2 — e ip 
perturbation theory, the coefficient of the eigenvalue e ±,pn is uniquely fixed as 

(2»-2)! 

Cm = 



(2.10) 



(n-iy.nl 

As we shall explain below graphically, e pn is the eigenvalue of the maximal- shuffling op- 
erator on the spin chain lattice. From this argument, we conclude that the coefficients of the 
maximal- shuffling operators are fixed by the assumption of the integrability and the represen- 
tation theory of off- shell symmetry superalgebra. In this work, we shall explicitly compute 
these coefficients at six loops and compare with (12.10I ). It constitutes a nontrivial check of the 
quantum integrability of 5\£ = 6 Chern-Simons theory. 

To derive the spectrum, following [35 J, we use the lattice-momentum eigenstates of elemen- 
tary magnon. For the A-spin magnon propagating on the odd-site chain, we have 

L 



\P\ 



L-\ 

£=0 



,..(2£+l) A2 ...) 



and 



\ph 



£=l 



Hp 



..(2£) Bl ...). (2.11) 



Here, | . . . (2£+ 1)a 2 • • •) refers that we put A2 at the 2£+Vs site while we put A\ for the remain- 
ing odd-sites, and similarly for the even-site chain. Hence, for an elementary magnon state, 
we may consider two kinds of states \p)a ® |0)fi and \0)a ® \p)b- Below we shall focus on the 
odd-site chain as the odd-site and even-site chains behave independently for the above subset of 
magnon states we are interested in. 

The corresponding integrable Hamiltonian at each order is well known for the above set of 
states 3. Since the structure of the even chain is identical to the odd chain, it is sufficient to focus 



3 All order generalization to the full spin chain that takes account of 4 — 4 interactions would be extremely 
interesting. 
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on the odd- site chain only. The zeroth order spin chain Hamiltonian counts the classical scaling 
dimension of the spins: 



1 



(2.12) 



1=0 



The two loop part of the Hamiltonian is given by [|2Tll22l 



L-l 



#2 = 4£ 



where 



%2 



[I -P : 



2i+l,2i+3. 



(2.13) 



£=o 



and P is the permutation operator defined by Wi = 8f S£. We shall take the infinite volume limit, 
L — > oo, and do not consider wrapping interactions. The corresponding 4-loop Hamiltonian can 
be identified as 



L-l L-l j 

#4 = e 4 , 2 £ <Q)2f 2 + e 4A £ Oft where Oft = — 4(1 - ¥ 2e+lj2e+3 ) - (I - P^+i^+s) 

^=0 £=0 



16 L 



This parallels the analysis of the 5\£ = 4 super Yang-Mills theory 
nian can also be identified as 



The 6-loop Hamilto- 



H 6 = e 6 ,2 £ + e 6A £ Oft + e 6 , 6 £ [( 1 - K 6 ) 0|f 6 + K 6 Oft 



L-l 

=0 



L-l 

E 

£=0 



(2.14) 



Here, K6 is an arbitrary coefficient. We thus see that a new feature arises beginning at six loop 
order. Up to four loops, candidate spin chain operators consistent with the quantum integrability 
are uniquely fixed. At 6-loops, operators consistent with the integrability are no longer unique: 
there are two commuting, mutually independent operators 



,2£ 
%.6 



,2£ 
%,6 



1 

64 



64 



1*21+1,2^+5 ^2£+3,2£+7 ~ P^+1,21+7 ^21+3,21+5 
+ 4P2£+l,2£+5 - 14P2£+l,2£+3 + 101 
; I - ^2£+ 1 ,2£+7 ) — 6(1 — F 2l+ j o e+5 ) + 1 5 ( I - P 2 ^+ ,2|+3 ) 



(2.15) 



Notice that these operator contents are identifiable with the Inozemtsev spin chain system 11381 
[8). Acting on the momentum eigenstate \p)a in (12.1 II) . these operators are diagonalized with 
the same eigenvalue 



£<2n 

t=0 



£=0 



\p)a = sin2 "f \p)a- 



(2.16) 
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So, the operator content of the six-loop Hamiltonian is not determinable uniquely by probing 
a single magnon dispersion relation. On the other hand, the operators (0>6,6 and ©6,6 are dis- 
tinguishable by acting on a state containing two or more magnon excitations. Note that both 
of them include the maximal-shuffling operators that produce the eigenvalues e ±inp . From the 
purely integrability point of view of the long-ranged Heisenberg spin chain, both operators are 
allowed and fit perfectly. Below, we shall approach this issue by computing relevant Feynman 
diagrams explicitly. We then determine the coefficient K6 and, from it, the operator contents of 
the spin chain Hamiltonian. 

3 Recursive Method for Maximal-Ranged Interactions 

We begin in this section with a recursive method, first put forward by Gross, Mikhailov and 
Roiban J3]|, of extracting maximal-shuffling terms among maximal-ranged interactions. We 
first redo the computation in 9£ = 4 super Yang-Mills theory and then repeat the computation 
in 5\£ = 6 Chern-Simons theory. 

Of maximal-ranged interactions, we are particularly interested in the coefficients Cjn of the 
maximal shuffling operators — they will provide a direct test of the integrability. In 9£ = 4 
super Yang-Mills theory, Feynman diagrams contributing to maximal-shuffling term are easily 
identifiable. As will be briefly reviewed below, at each order in perturbation theory, there is only 
one type diagram contributing to maximal-shuffling term. We refer to it as homogeneous dia- 
gram. For example, at 3-loop order, this maximal-shuffling term is responsible for the operator 
Og 6- There are also non-maximal-shuffling terms generated from maximal-ranged diagrams. 

For fA£ = 6 Chern-Simons theory, despite different interaction structures, we find a strik- 
ingly similar pattern repeated. Among the maximal-ranged interactions, maximal- shuffling 
term arises from homogeneous diagram. It turns out this diagram gives rise to the operator 
06,6, as in the situation of the 2\£ = 6 super Yang-Mills theory. We shall further confirm the 
coefficients C2 n of the maximal shuffling terms. Among the maximal-ranged interactions, there 
are also inhomogeneous diagrams that give rise to nonmaximal- shuffled terms. At six-loop or- 
der, for instance, these diagrams are responsible for the operator K6 [0*6,6 — 0*6,6 ]• However, 
by a direct computation, we shall find the coefficient K6 vanishes identically. This is interesting 
since, a priori, this operator could be generated given that interactions in the 5\£ = 6 Chern- 
Simons theory are different from those in the i\£ = 4 super Yang-Mills theory, this operator 
could be generated. 



9 



3.1 9t = 4 super Yang-Mills theory 

We first rederive the all-loop, maximal ranged interactions in 9£ = 4 super Yang-Mills theory 
obtained by Gross, Mikhailov and Roiban [3], emphasizing aspects directly relevant for similar 
considerations in ts[ = 6 Chern-Simons theory. 



Figure 1 : Homogeneous diagram in maximal-ranged interactions inO\[ =4 super Yang-Mills theory. It 
generate maximal-shuffling term in the dilatation operator. 

Consider in the 5\£ = 4 super- Yang-Mills theory renormalization of composite operators. 
Introduce in the theory the dimensional regularization parameter £, d = 4 — 2e of ultraviolet 
divergences. To avoid infrared divergences, a set of external momenta q are also injected to the 
operator. Denote by /bare(£ 5 ^) a collection of regularized Feynman diagrams that contribute to 
the renormalization of composite operators of the type (12.11) . Omitting wave function renormal- 
ization part proportional to the identity part, multiplicative renormalizability of the operators 
asserts that 

r 1 r 1 d? / \ i 

W?) = ex P [ Yz J J \ H4 W ~ 1 J J /bare ( e ' q ) (3 ' 1 } 

ought to be finite in the limit where £ goes to zero. The generator H^{X) is the quantum dilatation 
operator generating renormalzation group transformation and its eigenvalue corresponds to the 
magnon dispersion relation in the spin chain interpretation. In perturbation theory, 

oo 

We,9) = l + £/£L(e,9)- ( 3 - 2 ) 

i=\ 

From the £-independence of I K n(k), we get the relation 

/ ^(H 4 (kt)-l) =-lim 2£ln/ bare (£,4), (3.3) 

JO t V / e->0 
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for the operators which do not include the identity part. In planar 2\£ = 4 super Yang-Mills 
theory, the maximal-ranged interactions are generated only by quartic scalar interactions. Ho- 
mogeneous diagrams among them are depicted in Fig. [Q To avoid infrared divergences, we 
inject a finite momentum q to the X field from the right in Fig. [Q It suffices to keep zero mo- 
mentum for all other Z fields. The £-\oop contribution to the maximal shuffling diagrams can 
be evaluated recursively from the (£ — l)-loop contribution J3): 

i& 1 [r/1 „\ i£f l 



At) 

bare 



4% 



1 



[r(i 



T(l+ie) 



(3.4) 



efiu r(2-(/+l)e)n^(l-je) ' 

where a denotes X (e ip + e~ ip — 2) /(16ft 2 ) and p denotes the lattice momentum such that e ±,p 
generates shift one lattice site to the left or right in the spin chain. From the integrability, we 
expect the Hamiltonian to be 



H 4 (l,p) = Vl^4a = \ 1 



A% 2 



( e i P + e -i P _2). 



(3.5) 



One can check finiteness of the renormalized diagrams I ien (q) in (13.11) order by order in 
X. For instance, we checked this explicitly to the order 0(k 5 ) using the Mathematica and 
found that the renormalized diagrams to this order are indeed finite, canceling all £ " (n > 
0) singular terms. As was done in 0, we now show asymptotically that the singular terms 
of the regularized amplitude 7bare(£) are canceled by the expected Hamiltonian (13.51) in the 
renormalization factor 



exp 



_L /' 1 ^( v /T34^7_r 



2eJ t 



exp 



4a- 1+ln 



(3.6) 



in the limit £ goes to zero. To show this, we take £ — > while holding x : 
all-loop contribution by the Euler-McLaughlin formula: 



£ finite and sum the 



4are(e,<?) 



1 



{i) (e a) 

barel fc 'yJ 



Ve Jo 



1 



— / dxf4(x,q)exp - ( x(lna — lnx +2) + (1 — x) ln(l — x) 



We have relegated all sub-leading remainder to fa{x,q): 



U{x,q) 



J_ e -W)x I ^ 



47tVr(i+jc) 



m-x) 



.(3.7) 



(3.8) 



where \\f(z) is the poly-gamma function. The integral can be evaluated by the saddle-point 
approximation. At the saddle-point 



x = ^(i-VT^4a). 



(3.9) 
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the integral is evaluated as 

hanfaq) = ex P 



4a 



(3.10) 



1 + V1 

The £ 1 pole term is precisely inverse of the renormalization factor (13.61) dictated by the inte- 
grability. As it should be, the renormalization factor is independent of the infrared regularizing 
momentum q - the dependence resides in the finite remainder function R^{q). 



Figure 2: Inhomogeneous diagram in maximal-ranged interactions in 9£ =4 super Yang-Mills theory. 
It generate nonmaximal-shuffling term in the dilatation operator. The circles mark presence of vertices 
none of whose legs are connected to the operator. 

Among the maximal-ranged interactions, there are also inhomogeneous diagrams. At four- 
loop order, they were studied in 113411 . These diagrams are distinguished from homogeneous 
diagrams that some of the scalar quartic vertices do not connect to the operator at all. They are 
responsible for generating nonmaximal-shuffling terms in the dilatation operator. One can con- 
vince that these diagrams proliferate rapidly at higher orders in perturbation theory and, even 
worse, do not show any recursive pattern. Therefore, their contribution needs to be computed 
individually. In section 6, adopting two-point function method, we will find that the inhomoge- 
neous diagrams can be computed without ambiguity. 



3.2 9i = 6 Chern-Simons theory 

We now extend the recursive method to 3\£ = 6 Chern-Simons theory. One easily see that rele- 
vant Feynman diagrams are classifiable again into homogeneous and inhomogeneous diagrams. 
The homogeneous diagrams are planar irreducible diagrams all of whose interaction vertices 
are connected to the operator by two internal lines. The inhomogeneous diagrams: planar ir- 
reducible diagrams some of whose interaction vertices are connected to the operator by one 
internal line or not connected to the operator at all. 

Here, we first consider the homogeneous diagrams as they are easier to evaluate. The ho- 
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Figure 3: 'homogeneous' maximal shuffling diagrams 



mogeneous diagrams include all maximal shuffling terms, let us first state the expected scaling 
of operators via the anomalous dimension. 

Here, as explained in detail in G2ll25ll , we adopt dimension reduction 2(0 = 3 — £ for regu- 
larizing ultraviolet divergences in Feynman diagrams. To avoid infrared divergences, we again 
inject momenta q to the composite operator. Multiplicative renormalizability of composite op- 
erator asserts that 



W<?) = ex P 



1 



dt 



2eJo t 



H 3 (l 2 t) 



-fbare(£,<?) 



(3.11) 



is independent of the regulator e. Recall that, in 3\£ = 6 Chern-Simons theory, ultraviolet diver- 
gences arise only at even loops. So, in perturbation theory, 



We,<?) = l + £ J b2( e >?) 



(3.12) 



l=\ 



where q) denotes the 2£-loop regularized diagrams of order 0(k 2( ) . From £-independence 

of Iten(q), we get the relation 



/ ^ (m(\ 2 t) - i) = - lim 2£ In /bare(e,^) 
Jo t \ 2/ e^O 



(3.13) 



for the operators which do not include the identity part. 

As we explained in the previous section, from the integrability, we expect that the Hamilto- 
nian for the maximal shuffling term is given by 

H 3 (X 2 ) = ^ Vl-4b. (3.14) 

Here, b = X 2 (e ip + e~' p — 2) and e ±ip should now be interpreted as a left/right shift opera- 
tor by one lattice spacing on even- or odd-sites of the alternating spin chain. Therefore, the 
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integrability asserts that the renormalization factor is given by 



exp 



2eJo t y 2 2' 



exp 



r 1 / / ~ 2 

2e V 1 + VI — 4fo 



(3.15) 



To check if ultraviolet divergence of the bare diagrams /bare is inverse of (13.151) . we 
now evaluate the homogeneous diagrams. At elementary 2^ + 2-loops, from Fig. 2, we see that 
the diagram can be evaluated using £ many skeleton propagators defined by 1 -bubble diagram: 

one skeleton propagator defined by 2-bubble diagram: 

rm gg . r r 3 ((D-l)r(3-2(D) 

L3W= (Fp WhCre C3= (4^)^(3(0-3) (3 - 17) 

and using recursively the skeleton 1-loop integral: 

d 2(0 k 1 



G(a,b) := (4n) a (p 2 ) a+b - & f ■ 



(2%) 2(S > (k 2 ) a ((k + p) 2 ) b 
T(a + b- co) r(oo - a)r(co - b) 



(3.18) 



r(a)r(&) T(2a-a-b) 

Denote by the coefficient that skeleton £-loop contributes to the permutation P in the maximal- 
shuffling term = 21 — 4P in the dilatation operator. Label n-th loop momenta by k n and inject 
an external momentum q at the last skeleton vertex to regulate the loop integrals. At other 
vertices, there is no need to inject external momenta. Taking vertex and symmetry factors into 
account, the skeleton 0-loop (which is actually elementary 2-loop) contribution reads 

*,(*) = (-(2k) 2 ) -(-4)^^. (3.19) 

The skeleton 1-loop (which is actually elementary 4-loop) contribution reads 

d 2a k. 



2)2-<B 

= (-(2^) 2 ) 2 (-4) 2 G(4 - 200, 2 - «,) ■ 0-20) 
The skeleton 2-loop (which is actually elementary 6-loop) contribution reads 

Mi) = (-W 2 )(-4)/^ife)]g ((fe _V» 

= (4nf G(4 - 200, 2 - oo)G(7 -400,2- oo) ( ^ )2co ^ 3 3(3 _ 2a)) ■ 0-21) 
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Recursive pattern is evident. The skeleton £-th loop (note that this is actually elementary 2£ + 2- 
loop) contribution reads 



Mo) 



fjG(l+n(3-2(o),2-oo) 



n=l 



(47i) to (^ 2 )^ +1 )( 3 - 2(0 ) ' 



(3.22) 



Multiplying b e+l to Rf , we find the regularized 2£-loop amplitude has the expression 



r (2Q 

bare 



4tiV £ i [r(i-|)] 2m r(i+£ £ ) r b 



An 



(3.23) 



Using Mathematica, we checked up to 10-loop orders that the renormalized diagram I ien (q) in 
(13.1 II) with the H3 dictated by the integrability is indeed finite in the limit where £ goes to zero. 

As in the 5\£ = 4 super Yang-Mills theory, we can estimate the asymptotic behavior of the 
bare diagram 7bare(£, <?) m the limit £ — > °°, £ — > while holding x := £z constant. By the Euler- 
McLaughlin formula, we have 



1 



y/e Jo 



1 



2£ 



lb are (&,q)^^ dxf 3 (x,q)exp — [2x(ln b -In 2x +2) + (1 - 2x) ln(l - 2x) 



,(3.24) 



where fj (x) denotes a sub-dominant remainder 

1 



2x 



47tV r(i+x) 

7 ) fa^x) 



(3.25) 



We evaluate the integral by the saddle-point approximation. At the saddle point: 



1 . 

x Q = - (1 



1-46) 



the integral is given by 

W(e,<?) =exp 



1 

2i 



l-4£-l+ln 



+^ 3 (^) 



(3.26) 



(3.27) 



1 + Vl-46 

We see that the £ _1 pole term is precisely the inverse of the renormalization factor in (13.151) . It 
is remarkable that this all-loop agreement between the homogeneous diagrams and the integra- 
bility is closely parallel to the situation in the 5\£ = 4 super Yang-Mills theory. Because of this, 
in the next section, we adopt the two-point function method for deriving quantum dilatation 
operator. 

We also need to take account of inhomogeneous diagrams. It is easy to see that this class 
of diagrams does not include the maximal- shuffling terms. Nevertheless, the interaction range 
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Figure 4: 'inhomogeneous' maximal ranged diagrams 



is still maximal, viz. maximal-ranged, at a given order in perturbation theory, as illustrated in 
Fig.Sl In so far as one just focuses on maximal-shuffling part the spectrum, the homogeneous 
diagrams are sufficient. If one would like to identify operator contents of the spin chain Hamil- 
tonian, however, it is indispensable and crucial to take account of the inhomogeneous diagrams. 
For instance, at 6-loop order, the inhomogeneous diagrams are responsible for the coefficient 
k 6 of the operator J 6 n 6 = [6 6j6 - 6;6 ] . 

One would like to see if all-loop contribution of the inhomogeneous diagrams is also ob- 
tainable from the recursive method, much as for the homogeneous diagrams. Here, as in fA^ =4 
super Yang-Mills theory, the relevant diagrams proliferate rapidly at each higher order in per- 
turbation theory and do not exhibit recursive pattern in any obvious way. 



4 Anomalous dimension matrix 



From now, as in the previous works ||22[ 1251 [351 . we shall extract the anomalous dimension 
matrix of the single-trace operators O of the type (12.11) from two-point correlation functions: 

( : O(x)::O(0) : ) e = _f£ (4.1) 

From the dual Type IIA string theory viewpoint, this method amounts to deriving time-evolution 
Hamiltonian of a single non-interacting string propagating in AdS4 x CP 3 spacetime. We use 
the dimensional reduction to regularize ultraviolet divergences. The two-point correlation func- 
tions are related to Feynman loop diagrams A21 by 

(: 0(jc) : : O(0) :) e = (7 e ) 2L g-fcM* 2 A 2 (*0) 

= (7 e ) 2L exp [ ln(l +A 2 X 2 +A 4 l 4 +A 6 l 6 •••)], (4.2) 
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where 7 e denote the Euclidean scalar propagator in the position space 



(2tc) 2co p 2 



4K a (x 2 ) 



(4.3) 



Here we consider all Feynman diagrams, connected or not, contributing to a given order 
of X 2 ^. As the definition of the anomalous dimension matrix takes the logarithm, it suffices 
to compute connected diagrams only. The anomalous dimension matrix is then extractable as 
coefficient of ln(;c 2 ) in the exponent. Up to 6-loop orders, the spin chain Hamiltonians // 2 £ 
classified in section 2 are given by 



H2 = — lim £ A 2 



H4 = - lim 2e 



H 6 



• lim 3£ 



1 1 1 

A 6 --(A 2 A 4 +A 4 A2) + -A i 2 



(4.4) 



Note the extra factor n multiplied. It arises combinatorially from extracting coefficients of In x 2 
from the 2£-loop contribution (lnA) 2 ^, 



(lnA) 2l = [c2ge~ l +0(e°)] (x 2 ) ee = \c u e 1 + 0(e°)] [l +£e hue 2 + 0(e : 
At 2£-loop order, we have 



Hip = —£ lim £ (In A)™ . 

e^O 1 



(4.5) 



(4.6) 



Summing over all loops, 
and we have 



and A(£) = l + £^A 2/ (e) 
l=\ 



(4.7) 



-lim £ lnA(e) 



lim 2£ In 7 bare (£) 

e-j'O 



(4.8) 



Thus, starting from a different prescription, we have arrived at the same definition of the dilata- 
tion operator as (13.131) . Here, the last equality followed from the fact that the A(e) is defined in 
terms of the two-point correlation function ( :0 (x) : :0 (0) : ), so is square of the 7bare(£) (equiv- 
alently, the lnA(e) is twice of the ln/bare(£) connected diagrams) we considered in subsections 
3.2 and 3.3. 
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Multiplicative renormalizability of the composite operators asserts that all £~ n , (n > 2) sin- 
gularities must cancel each other. We shall use such cancellations as a checkpoint of our com- 
putations. In A2, the contribution starts from e power and H2 is automatically finite. In A4, 
leading singularity starts from £ 2 power. The coefficient of this singularity in (lnA)4 should 
vanish and the coefficient of order £~ ! leads to H4. In A$, there are in general singularities of 
£~ 3 and £ 2 powers. The coefficients of them in (InA)^ vanish and the coefficient of order £ _1 
leads to H^. As discussed in Section 2, beginning at 6-loop orders, several independent spin 
chain operators can appear in H. The cancelation of higher-order singularities must then take 
place for the coefficients of each independent operators. At 6-loops, this will provide an addi- 
tional stringent consistency check of our earlier 4-loop results ll35l on whether computation of 
A4 was correct and identification of spin chain operators was complete. 



Figure 5: 6-loop homogeneous diagram involving maximal shuffling terms 



Figure 6: adjoint of the above homogeneous diagram 



5 Operator Contents of 6-Loop Dilatation Operator 

We asserted earlier that, beginning at 6-loop order, the dilatation operator becomes complicated 
because there arise two types of maximally ranged interactions. In this section, we shall explain 
this in detail by analyzing operator contents of the 6-loop dilatation operator, focusing on the 
maximally ranged diagrams involving seven lattice sites. 
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Figure 7: 6 loop inhomogeneous diagram whose interactions are maximal-ranged 



Figure 8: adjoint of the above inhomogeneous diagram 



The first type of spin chain operator arises from the homogeneous diagrams and they give 
rise to the maximal shuffling terms. The diagrams in Figs. \5\ and [6] belong to this category. The 
sextet interaction vertices are combined as in Figs. [5] and [6] The second type arises from the 
inhomogeneous diagrams. They are still maximally ranged but do not generate the maximal 
shuffling terms. Figs. [7] and [8] depict the inhomogeneous maximal ranged diagram at 6-loop 
order. 

Their operator contents can be constructed from the basic building structure arising from 
the 2-loop sextet scalar vertices: 

Of 23 = 21 — IQ + i^ + 2 — Kj> +2 j + 3 + ^£+l,£+3^i+l,£+2 

+ 2P i+M+3 ]EQ + 2,*+3-4P/+ M+ 3, (5.1) 

where K^' = § a t,h a ' b ' is the 4—4 co variant contraction operator. Here £ could be even or odd 
but, as mentioned before, we focus on the odd chain where £ is even. Acting this operator on the 
set of states in (12.21 ), all the terms involving the contraction operators drop out and the operator 
is reduced to 

05 23 =2I-4P/ +M+3 . (5.2) 
Below we shall present this reduced part only and omit all the contraction type operators. The 
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operator structure of the diagram in Fig. [5] is then 

®123 O345 °567 = 2 3 (I - 2P* +M+3 ) (I - 2Pf +3 ,f +5 ) (I - 2F e+5j£+1 ) 
Similarly, the operator structure of Fig. [6] becomes 

O e 567 O e 345 Of 23 = 2 3 (I — 2¥ e+5 j +1 ) (I - 2¥ e+3J+5 ) (I - 2F i+u+3 ) 



(5.3) 



(5.4) 



The numerical parts of these Feynman diagrams are the same. So, adding the two contributions, 
we obtain the spin chain operator 



>?~7 



il3 0345 0«7 + 0»7 °345 ^ 23 . 



One can show that this operator is related to the 6 by 

£ Ot 7 = 64 £ \640{ 6 - S0O{ 4 + 220{ 2 - - 1 



(5.5) 



(5.6) 



The details of the derivation is relegated to appendix B. We focus on the maximally ranged 
interaction part only, and omit interactions of lower range given by e6,2 ©2,2 and ^6,4 04,4. Thus, 
we conclude that the operator content of the homogeneous diagram is given by 



®L 7 = 4 6 ol 6 +--- 

where the ellipses denote the interactions of lower range. 

The operator structure of the diagram in Fig.|7]is identified as 

O e 567 0{ 23 0/345 = 2 3 (I — 2P* +5 ,/+ 7 ) (I - 2P/ +M+3 ) (I - 2F e+3 , e+5 ) 
and that of the diagram in Fig. [8] as 



(5.7) 



345 Of 23 ©567 = 2" (I - 2P £+3) £+5) (I " 2P* +M+3 ) (I - 2P /+5 ,£ +7 ) 



(5.8) 



(5.9) 



Notice that the order the basic unit operators are multiplied is different from that in the ho- 
mogeneous diagrams. Again, the numerical parts of these Feynman diagrams are the same. 
Combining them together, we obtain 



»L 7 



In appendix B, we show that 



»L 7 



^OLsO^O^ Of 23 04,7. 



C 7 -64P /+M+7 + 32I 

C? + 64 [ 64 ©6,6 " 96 ©4,4 + 36 ©2,2 " 2 1 
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(5.10) 



(5.11) 



Hence, we have 



>L 7 = 4 6 [- 



+ 



(5.12) 



where ellipses denotes the interactions of shorter range given by e6,2 ®2,2 and ^6,4 O4 4. 

For completeness, we review the operator structures arising in the 4-loop maximally ranged 
interactions. For the 4-loop order, only diagrams of homogeneous type are present. These Feyn- 
man diagrams are depicted in Fig. |9] and Fig. [Kjl The numerical parts of the 4-loop amplitudes 
in Figs. [9] and [TO] are the same. Combining them, we have the spin chain operator 



'1-5 
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>123©345 +®345 ®123 



-I — £+3 — P^ + 3 £ + 5+F£ + i £ + jF£ + 2 £+5 +Fi + s t+sFt+i £ + 3 



This can be rewritten as 



(5.13) 



i6P t+u+5 -m 



(5.14) 



where we again omit the contraction type operators. Here, we used the identity 

ei ai a 3 a 5 £lb lb3 b 5 = I + Pl3 P35 + P35 Pl3 " Pl3 " P35 " Pl5 = , 

taking account of the fact that at and b( are running over only {1,2}. 



(5.15) 



Figure 9: 4-loop homogeneous diagram involving maximal shuffling terms 



Figure 10: adjoint of the above 4-loop homogeneous diagram 
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6 Computation in the 6-Loop Dilatation Operator 



In this section, we shall compute maximal-ranged interactions in the 6-loop dilatation operator 
(14.41) and extract the coefficients ^6.6 an d Kg. As explained in the previous section, it suffices to 
compute the 6-loop maximally ranged operators, Og^ and Og g . 

The first contribution to the maximal-range interactions arise from the A\ part. The relevant 
2-loop amplitude in Fig.dUis organized as follows: 



A 2 = fl2 (e) {x 2 %f- 2& [ • • • + 0( 23 + O e M5 + ( 
In Ref. [35J, the numerical coefficient 02(2) was obtained as 

[r(o)-l)] 3 [r(3-2a))] 2 r(5a)-6) 

«2(£J = t 

1671 [r(3(0-3)] 2 r(6-4a)) 



^567 



+ • 



(6.1) 



1 

4i 



l + (l- V (I))£ + ^(ll7C 2 + 12 ¥ 2 (i)-24 ¥ (I)-72)£ 2 + 0(£ 3 ) 



(6.2) 



where the poly-gamma functions take the value V|/( j) = — C — 21n2 and V]/(l) = — C with C (= 
0.577215 • • •) being the Euler's constant. From this, we extract contribution of the A\ to the 
relevant maximal-ranged interactions as 



A 3 2 = (a 2 (e)) 3 (x 2 %f- 6a 
where we already omit all disconnected contributions 



+oi /7 +2GC 7 + . 



(6.3) 



Figure 1 1 : 2-loop diagram involving three sites 

The (A2A4+A4A2) part is another contribution to the maximal ranged interactions. The 
relevant 4-loop contribution to the maximal-ranged interactions has the structure 



A4 = (X 71) 



+ 04(e) [OU + 0^7 ] + 4(e) Q{ 23 O e 561 + 



(6.4) 



This along with (|6.1I) gives rise to six loop maximal-ranged interactions through (A2A4 + 
A 4 A 2 ). We find 



A 2 A4+A 4 A 2 = (;c 2 7t) 9 - 6co 



+ 2a 2 (£) 04(e) [ 0^ 7 + Of „ 7 ] + 4(e) Q^ 7 + ■ 



(6.5) 
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Lastly, as discussed before, the relevant maximal-range interactions in are given by 



A 6 = (x 2 %) 9 - 6a 



(6.6) 



Using (14.41) . we found that the coefficients of the maximal-ranged operators in (lnA)g has the 
structure 



(lnA)e 



+ 



1 , 

a 6 (e) - a 2 (e)a 4 (e) + ^a 2 (e) 
1 , 

a 6 (£) - a2(e)a4(e) + -ra^e) 
6 



^1~7 



»W7 + ' 



(6.7) 



The 4-loop coefficient 04(e) was determined in 1135 



a 4 (e) =7i(e) 



:r(o>-l)]- 

4 4 7l 2 



r(2-co)r(3-2a)) 
J(2©-2)r(3(D-3) 



r(9co-n) 
r(n-8co) 



(6.8) 



The integral J\ (e) is introduced by the following skeleton 2-loop integral K5: 

d 2a k d 2a l 



(2tc) 2co (2tc) 2co (fc 2 )^(/ 2 ) 



W'2 



1 



(6.9) 



~[(i-p) 2 P[(k- P ) 2 p[(k-i)T 5 - 

The corresponding skeleton 2-loop integral is depicted in Fig.Q/l] As we evaluate the two-point 
correlation functions in the x space, we first evaluate the skeleton graph with nonvanishing 
momentum injected from the left to the right and then Fourier transform it back to the x-space. 
See Ref. [35] for the details of the method. The integral J\ (e) evaluated from the skeleton 2-loop 





Figure 12: The left is a skeleton 2-loop diagram used to evaluate A4. The skeleton 4-loop diagram on 
the right is for A^. We inject a momentum p following the arrows in each diagram. 

integral ^5: 



7i(e) =^ 5 (2 -co, 3 -2(0, 2 -(0,3 -2co, 1 



(6.10) 



23 



This integral is already nontrivial and does not allow analytic evaluation. Here we evaluated it 
numerically by using a Mathematica package. We took the strategy of expanding the integral 
7i (e) in powers of £ and determine term by term. The result is 



/i(e) 



a 



l+a^ 1 £ + a) 1 £ 2 + --- 



371 £ 
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1 + ( 1 3 + V|/( \ ) ) £ + ( 488 - y n 2 + 52\|/( \ ) + 2y 2 ( 



+ • 



(6.11) 



We computed the leading coefficient cc^ 1 numerically and matched to the above rational value. 
This is the value required for the multiplicative renormalizability of composite operator at four 
loops. The second coefficient is the one relevant for the 4-loop amplitude 0Cy . Here, we checked 
this against the requirement of the multiplicative renormalizability at six loops. We computed 
the last coefficient ocj numerically up to 30 significant digits. We match this number as a linear 
combination of 1 and transcendental numbers 7C 2 , ^(3), y (j) with the requirement that the 
coefficients are rational numbers. This is highly nontrivial to check. However, utilizing the 
remarkable algorithm PSLQ discovered recently by Ferguson, Baily and Arno, we determined 
the combination in (16.1 II) as the unique solution. The details are explained in appendix D. Using 
this, we found that 614(E) has the series expansion 



04(e) 



32e 2 



1+ 3-2 ¥ (i) 8+ 



2 + ^7r 2 + 2 ¥ 2 (i)-6 ¥ (I) 



£ 2 + . 



(6.12) 



We now turn to the 6-loop coefficients a§{€) and d^(e). Consider first the coefficients arising 
from the diagrams in Figs. [5] and [6l One can easily verify that these two coefficients are the 
same. As explained in appendix C, we find 



a 6 (e) 



J 4 {e) [r(o)-l)] 7 [r(2-o))] 4 [r(3-2o))] 2 r(13(o-16) 



4 6 7C 3 



[r(2(0-2)] 4 [r(3(0-3)] 2 r(16-12(o) 



(6.13) 



The Feynman integral can be reduced analytically to the function /4(E). Consider the skeleton 
4-loop integral defined by 



£ 8 (wi---w 8 ) = (4tc) 4co ( p 2 )- 4ro +^f=i^ j 



d 2a h d 2a k 2 d 2a k 3 d 2(0 k 4 
(27t) 2co (27T) 2m (27t) 2co (2tt) 2co 



1 



(k\)^ (k 2 2 yi (k 2 )^ (k 2 y^ [(k l+ k 2 - P ) 2 p 

1 

i(k 3 +k 4 -p) 2 p \{h-h) 2 r [(* 2 -*4) 2 P 



(6.14) 
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Then the function /4(E) is given by 

y 4 (e) =# 8 (2-(0,3-2(D,3-2(0,2-(0,2-(D,2-(0, 1,1). 



(6.15) 



By using the PSLQ algorithm again, we determined that the function /4(e) has the expansion in 
powers of £ as 



16 



135tt 2 £ 2 



/ 14 \ 
l + (28 + - + 2 ¥ (I)j£ 



+ 



c ^ 11 2 27 2 868 n . „ 2/1< 

522 ~ 25 -9 - y 71 +75-V(5) + 2^ (i 



(6.16) 



The coefficients of the first two terms fit precisely to the values required for the multiplicative 
renormalizability. As presented in appendix D, we checked numerically the first coefficient up 
to 14 significant digits and the second up to 9 significant digits. From the expansion (16.161) . we 
found that a^E) is given by 



a6(£) 



1 



384£ 3 



80 11 



1+3- 2- ¥ (i) £+ - + ^ + 4 TO-16¥G) e 



Adding all pertinent contributions to Og 6 , we get 

1 ^ 

a 6 (e) - a 2 (e)a 4 (e) + -a|(e) 



96£ 



+ 0(e°). 



. (6.17) 



(6.18) 



Contributions to the £ 3 and £ 2 terms cancel out, thus satisying the consistency of multiplica- 
tive renormalizability. From (14.41) and (15 .7fl , we finally determined the coefficient of Og 6 as 



e 6 ,6 = 128. 



(6.19) 



This is the value that agrees exactly with the prediction of the integrability. 

Similarly, the coefficient d^it) in (16.61 ) arises from the diagrams in Figs. [7] and [8j As ex- 
plained in appendix C, we find 

/ 5 (e) [r(co-l)] 6 [r(2-co)] 2 [r(3-2a))] 3 r(13(0-16) 



4 6 rc 3 [r(2co-2)] 2 [r(3co-3))] 3 r(16-12(o) 
where the function /5(E) is the skeleton 4-loop integral: 

7 5 (e) =K%(2- (0,2 -00, 3-2oo, 3-2co, 3-2oo, 1,1,1). 

Again, utilizing the PSLQ algorithm, we computed the series expansion of 75 (e) as 



(6.20) 



(6.21) 



/s(e) 



16 



/ 14 9 \ 
l+(l9 + - + - + ¥ (l)+x|/(I)j8 



+ 



135e 
177587 61 



733 



im2 



225 



6 



+ ■ 



.(6.22) 
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We confirmed that the first two coefficients are the values required by the multiplicative renor- 
malizability. See appendix D. We then obtained the coefficient a^(e) in £ expansion as 



a 6 (E) 



1 



192£ 3 



11 



l + ( ^-3 V U))e+( f ^+| V 2 (i)-y ¥ (i) )e> + 
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(6.23) 



From this, we find the pertinent contribution to Og 6 as 

a 6 (e) - a 2 (£)a 4 (e) + ^a\{z) = • ^ + 0(e°) . (6.24) 

Again the coefficients of £ 3 and £ 2 are vanishing, which is the requirement of the renormal- 
izability. From (14.41) and (15 .120 , we finally find that 



K 6 = 0. 



(6.25) 



As such, we conclude the operator 0(,^, even though the integrability and the excitation sym- 
metry allow it to be present, does not contribute to the 6-loop dilatation operator. Curiously, 
this parallels exactly to the structure of the dilatation operator in 9t = 4 super Yang-Mills the- 
ory (HI — in both theories, the maximally-ranged interactions are governed by a single spin 
chain operator. 

Summarizing, we determined uniquely the operator contents and the recursive structure of 
the 6-loop dilatation operator given by 



fl 6 = 128jX 6 4- 



(6.26) 



where the ellipses denotes shorter-ranged terms we omitted. 



7 Remarks on AB J Theory 

Up to this point, we focused, among the fA£ = 6 Chern-Simons theory, on the ABJM theory. 
This theory has the gauge group U (N) x U (N), the Chern-Simons coefficients +k, —k and is 
invariant under parity. Notice that the parity is generalized to interchange the two gauge groups 
and hence all matter fields and their conjugates. 

We can extend straightforwardly our considerations to the ABJ theory, which has the gauge 
group U(M) xU(N), the Chern-Simons coefficients +k,-k but with N < M < N + k. This 
theory then breaks the generalized parity. As explained in [|25l and recapitulated in section 2, 
the integrability and the Yang-Baxter equations therein implied that the dilatation operator is in 
general parity non-conserving. Nevertheless, explicit computations in [|25l[35l showed that the 
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dilation operator sustains to be parity conserving. At 6-loop order and beyond, as there begins 
to arise two types of spin chain operators in the dilatation operator, it becomes interesting if 
they source a room for parity (non)conservation beginning at this order. 

From the Feynman diagrammatic s, however, it is evident that all the changes in 6-loop 
diagrams of the ABJ theory compared to those in the ABJM theory is rather trivial; one just 
replaces the 't Hooft coupling- squared X 2 of the ABJM theory by the product of two 't Hooft 
couplings XX of the ABJ theory. The nullification result K6 = in (|6.25l) of the ABJM theory 
was not sensitive to the 't Hooft coupling as its dependence is an overall weighting factor. In 
the ABJ theory, by repeating the color factor counting in the Feynman diagrams, we find that 
all the contributing diagrams have the common dependence on M and N. Therefore, K6 = in 
the ABJ theory as well, viz. there arises a unique spin chain operator to the six-loop dilatation 
operator. This then eliminates the possibility that, within maximal ranged interactions, parity 
non-conserving effect arises at six loops and beyond. 

On the other hand, we suspect that the effect of parity non-conservation may show up at 
planar limit in some of the higher conserved charges Q«, Q,„ (n > 3). With indications from 
recent results of [31 J, we also expect that parity non-conserving effects will also arise from 
non-planar corrections. Given that there is no known examples of parity non-conserving yet 
integrable system, further understanding on this issue would be very rewarding. 

8 Conclusions 

In this paper, we studied the dilatation operator of i\£ = 6 Chern-Simons theory, paying attention 
to two pertinent issues that begins to arise from six loops and beyond: magnon spectrum and 
operator contents of the dilatation operator viewed as a spin chain Hamiltonian. The integrabil- 
ity together with the excitation symmetry plus the 2-loop results led to the unique prediction for 
the maximal shuffling terms to all orders. 

We computed the maximal shuffling terms in the dilatation operator. We found that the 
coefficients of this term, which arise from so-called homogeneous diagrams, agree with the 
prediction of the integrability and the recursive property. We thus found that the situation of 
maximal shuffling interactions is exactly parallel to the situation in 3\£ = 4 super- Yang-Mills 
theory. Despite the parallel, we also argued that the 3\£ = 6 Chern-Simons theory is signifi- 
cantly different from the = 4 super Yang-Mills theory since there separately exist so-called 
inhomogeneous diagrams. These diagrams do not generate maximal shuffling terms but are still 
maximally-ranged. We showed that the inhomogeneous diagrams are not recursive and appears 
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to depend on the choice of infrared regularization. 

To handle these difficulties, we adopted ab initio method of extracting the dilatation oper- 
ator from two-point correlation functions of single-trace composite operator. In dual Type IIA 
string theory, this amounts to deriving time-evolution Hamiltonian of a single free string. The 
method was particularly convenient for extracting not only the spectrum but also the operator 
contents of the dilatation operator. By explicit computations, we determined the maximally- 
ranged interactions of the dilatation operator up to 6-loop order. We were able to compute 
pertinent Feynman diagrams by utilizing the remarkable integer-relation algorithm PSLQ. We 
found that, though details of Feynman diagrams contributing to the dilatation operator are very 
different, the operator contents and the recursive coefficients of the maximal-ranged interaction 
is identical to the dilatation operator in the 2\£ = 4 super- Yang-Mills theory. This result is of 
course consistent with the integrability, which we consider as a highly nontrivial test for the 
9i = 6 ABJM theory. 

We pointed out that extension of the results to the parity non-conserving 5\£ = 6 ABJ theory 
is rather trivial, at least for the maximally-ranged interactions. One just replaces the 't Hooft 
coupling squared X 2 of the fA£ = 6 ABJM theory by product of the two 't Hooft couplings Xk of 
the 3\£ = 6 ABJ theory. The resulting dilatation operator remains to be parity conserving. 

In the above test, we assume that the pseudo -momentum coincides with the lattice momen- 
tum in the operator space. This assumption can be further tested computing the next-maximal 
shuffling coefficient e^, which is predicted as ^6,4 = — 32/i2- 

Recently, the coefficient was computed in [|36l . Given that Ward identity was not verified 
to the order the computation is based on, the result of lf36l needs to be checked independently. In 
particular, there is no first-principle understanding of its origin. One promising approach of the 
independent check would be the following. In Ref. 11351 , it was noted that h(X) should also be 
present as a renormalization factor in the central charges of the off-shell \psu(2\2) © psu(2\2)\ tx 
IR 2,1 superalgebra. Any further understanding of this function h(X) such as renormalization of 
the central charges would be extremely interesting. 
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A Comparative calculation of inhomogeneous diagram 



Below, for comparison with two-point function method, we shall illustrate a comparative calcu- 
lation of a sample inhomogeneous diagram in operator-mixing method. Consider the inhomo- 
geneous diagram in Fig. 3. The diagram has three skeleton vertices but, by reflection symmetry 
of the diagram, there are only two independent skeleton vertices. To control potential infrared 
divergences, it is necessary to inject nonzero momentum to the operator judiciously diagram by 
diagram. Here, consider injecting q at the top or bottom skeleton vertex. This is what we did for 
the homogeneous diagrams and found to yield an infrared finite result. For the skeleton 2-loop, 
which is actually elementary 6-loop, we find the contribution 

1 



B 



bare 



= (-(27t) 2 ) 3 (-4) j C^G(l +6, 1)G(3 - © + £,£ 



1 



{Any 



ATI 



(4^)2co (^2)3e 

[r(i-f)] 7 r(i+3e) 

3!e3 (k-h)(k+h)(k-hnk-h: 



3e 1 



(A.l) 



Along with lower-range diagrams in (13.191 13.201) . this contributes to the 6-loop part, W(,, 

of ln/ bare (£): 

,-| bare 



W 6 (e) 



1 



1 



i) in 



(A.2) 



We relegate details of computation of the operator content to section 6 and focus on the coeffi- 
cient of Og n 6 . By multiplicative renormalizability, the leading singularity of must begin at 
order of 0(e.~ l ), viz. lim e ^o e m ^bare (e) must be finite. For this diagram, we found 



W 6 (e) 



1 , 25 + 9w(l) 2 . , s 
16 48 + 1 ' 



l 3 Oj, n 6 . 



(A3) 



Alternatively, we may opt to inject momentum q at the middle skeleton vertex. In this case, 
we find 
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'bare 



(-(2tt) 2 ) 3 (-4)^G(1 +£, l)G(2 + e-co, 1 + e) 
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2co („2\3e 



(47C) 2ro (q 
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3e 1 



[r(i-f)] 7 r(i + 3e) 
3£3 ( i-is) 2 (i-is)r(i-|: 



(A.4) 
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Hence the 6-loop part of the logarithm becomes 



(e) 



1 



bare 



•4 



O 

0.£- 3 + I e - 2 + O(£ _1 ) 
6 



3 /nun 



"6,6 

4 3 Oj, n 6 . 



(A.5) 



The Feynman diagram ought to be infrared finite once independent momenta are injected 
to every external vertices. Unfortunately, the loop integral in this case is too complicated and 
do not permit all-loop resummation. So, we see that in both options the problem stems from 
non-analyticity as some of the injected momenta are taken to zero. 



B Proof the operator identities d£6D and (|5.11[) 



To prove these two identities, we use the identity in (15.151) repeatedly. For (15.61) . we first show 
the following identity: 



Pl 7 P 35 - Pl5P37 = P13P35P57 + P57P35P13 " Pl3P57 " Pl5 - P37 + I • 



To show this, note that 



P 17 P 35 " Pl5P37 = (Pl5P57 + P57P15 " Pl5 " P57 + l) P35 



-Pl 5 (P 35 P 5 7 + P57P35 " P35 " P57 + I y 
-P35P13P57 + P57P35P13 " P57P35 " Pl5 + P35 + P15P57 
-P35P13P57 + P57P35P13 " P57P35 " Pl5 + P35 
+ (P13P35 + P35P13 " Pl3 " P 35 + l) P57 
P13P35P57 + P57P35P13 " P13P57 " Pl5 " P37 + I 
- (P35P57 + P57P35 " P35 " P57 " P37 + l) 
P13P35P57 + P57P35P13 " P13P57 " Pl5 " P37 + I , 



(B.l) 



(B.2) 



Here, we used (15.151) to replace Pn and P37 in the first line and the second P15 in the third line. 
Next, we turn to the operator 0^7: 

h230345©567 + 05670 3 450i23 
= (2Pi3-I)(2P 35 -I)(2P 5 7-l[) + (2P57-I)(2P35-I)(2Pi3-i; 
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= 8 (P13P35P57 + P57P35PB + P13P57 " Pl5 " P37 + I 

+4Pi 5 + 4P37-4P 35 -2I 
= 8 (P17P35 - P15P37) + 4P 15 + 4P 37 - 4P 35 - 21 



8 



- 640 6 , 6 + 80O 4 , 4 - 220 2 , 2 + ^-1 



(B.3) 



which is the proof of the identity (15.61) . 
Finally, we turn to the identity in (15.111 ): 

Oi^7 + Oi^7 = O123 (©345©567 + 5670 34 5 ) + (O345O567 + 5 670 3 45 ) 123 

= 2(I-2P 13 )(-8)(I-2P37) + (-8)(I-2P 37 )2(I-2P 13 ) 

= 32^I-2Pi 7 ). (B.4) 

C Derivation of ae (e) and de (e) 

The diagram in © involves four 1 -bubble and two 2-bubble propagators in total. The corre- 
sponding factors are (C 2 ) 4 (C3) 2 . This then leads to the skeleton 4-loop diagram in Fig. [121 We 
assign momentum and the w-factor to each internal line as follows; k\ and w\ = 2 — co for the 
upper-right line, p — k\—k 2 and w$ = 2 — CO for the middle-right line, fc 2 and w 2 = 3 — 2co for 
the lower-right line, £3 — k\ and wq = 1 for the upper-center line, k^ — & 2 and w% = 1 for the 
lower-center line, £3 and W3 = 3 — 2co for the upper- left line, p — k^—k^ and = 2 — co for the 
middle-left line, £4 and W4 = 2 — co for the lower-left line. One has then 

« 6 (8) (/ e ) W- 6co (4^55 • (C 2 ) 4 (C 3 ) 2 / {2n) l {p2y6 -n. *** > <P.l) 

where S = (— 4tt 2 ) 3 is the vertex and symmetry factors. 
For the Fourier transform, we used the formula, 



/ 



d 2 > 1 I>JC _ r(<D-a) 1 

(2tt) 2co (p 2 ) a 4 a 7T co r(a) (x 2 ) & - a ' 

This leads to the expression a(,(e) in (16. 131 ). 

The diagram in © involves two 1 -bubble and three 2-bubble propagators in total. This leads 
to the skeleton 4-loop diagram in Fig. [121 We assign momentum and their w-factor to each 
internal line as follows; k\ and w\ = 2 — co for the upper-right line, p — k\—k2 and W5 = 3 — 2co 
for the middle-right line, fc 2 and w 2 = 2 — CO for the lower-right line, £3 — k\ and wj = 1 for 
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the upper-center line, £4 — k 2 and w$ = 1 for the lower-center line, £3 and W3 = 3 — 2co for the 
upper-left line, p — kj — £4 and = 1 for the middle-left line, k<\ and W4 = 3 — 2(0 for the 
lower-left line. One has then 

1 r d 2a n 1 

36(e) (/ e ) V*) 9 - 6 " = S- ^ ./ 5 (e) • (C 2 ) 2 (C 3 ) 3 / (C3) 

This leads to the expression 05(e) in (|6.20l) . 



D Evaluation of the integrals J\ (e) , /4(e) and / 5 (e) 

In this appendix, we evaluate the higher order terms of the two-loop integral 7i (e) and two 
four-loop integrals, /4(E) and /5(e) in the section 6. It is difficult to find their analytic forms 
in a direct manner. We first numerically evaluate the integrals and then find the corresponding 
expressions consisting of transcendental numbers. There exists PSLQ algorithm[391 which is 
quite useful in finding such an analytic form from numerical data. For the numerical evaluation, 
we use the Mathematica packages-MB ||4~0l |4H and AMBRE [|42l . To find the corresponding 
analytic expression, we use the package PSLO.nb fl43l . which is a Mathematica implementation 
of the PSLQ algorithm. 

Since these methods are based on the Mellin-Barnes representation of Feynman integrals, 
we cast the integrals into the form: 

„/ s fdzi / , dz 2 r(-zl)r((D-W25-zi)r((0-wi+zi) mi . 

K 5 (w l ,W2,w 3 ,w 4 ,w 5 ) = 1 — 1 — — (D.l) 

J 2m J 2ni I (wi — zi) 

r(-Z 2 )r((0-W35 -Z 2 )r((0-W4+Z2) 



r(w 4 -z2) 

r(-(0 + VVi4 -Zl2)r(-(0 + W 2 35 +Zl + Z2)T(w 5 + Zn) 

r(w 2 )T(w 3 ) r(ws) r(w 2 35) r(2(o - ^235) r(2© - wu + zn) ' 

with W12- = w\ + W2 H and Z12... = z\ +zi H • Direct application of these packages yields 

the following result: 

1 



oc 

/x(e)= £5(2- a>,3 -2<o,2- (0,3 -2a>,l) = — 1 + a^e + a^e 2 + 0(e 



(D.2) 



where we have found cc^ 1 = — ^ and = 13+\|/(j) as reported in ll3~5l . The numerical value 



of ocj is obtained as 



4a), = 309.7165844821997330276736227359 .... (D.3) 
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By the PSLQ algorithm, we matched this numerical result as a linear combination of 1, n 2 , 
and \\f 2 ({) under the condition that their coefficients are rational numbers. This basis 
follows from requiring that e 4:4 is a rational number. This leads to the unique set of fractional 
coefficients: 



4< 



488 



17 



% 2 + 52\|/(i) + 2y 2 ( : 



(D.4) 



Note that the number of significant digits in this numerical value is 30 and the difference be- 
tween these two expression is less than 10~ 30 . We have checked the other numbers, a^ 1 and 
0$ to the same accuracy. 

The four-loop integrals Kg(wi, ■ ■■ ,w$) in (16.141) can be evaluated in a similar way. Its 
Mellin-Barnes representation is given as 



h(wi 



dzi 
2%i 



dzs N 
2ni D 



(D.5) 



where 



N 



and 



D 



Z246) T(W7 +Z235) r((0/2 + Z3g) T((0/2 — W48 — Zl) 

n=l 

) r((H)/2 ~W 6 - Zl) r(tO - W3467 - 2234) r((D - W 36 78 - Z256) 
) r(3(0/2-W334 67 7 8 -Zl223456)r(3(0/2-W 2 34678 -223457) 
) r(-2(0 + W 12 345678 +247)r(3d)/2-VV2345678 -Z23468) 
) r(2(0 — Wi2345678 —23478) r(—3d)/2 + W2345678 + 223478) 

) r(-(0 + W234678 + 22345678) T(-G) + W 346 78 +2123456) C°.6) 

= r(w 3 ) r(w 5 ) r(w 6 ) r(wj) r(w 8 ) r(to - w 46 8) r(w 5 - z 5 ) r(-z 38 ) 

<g> r(3©/2 - W34678 - Zl) r(2(D - W2345678 ~ 2234) 

(8) r(2(0-VV33466778 -2223456) r(-(0 + W234678 +223456) 

<g) r(5(D/2 - W12345678 - 247) T(-3(f)/2 + W12345678 + 23478) ■ (D.7) 



The 74(e) and 7s (e), respectively, in (16.161) and (16.211 ) are defined with help of Kg. By 
numerical computations, we find that 7 4 (e) has the expansion, 

0.012008436579832... 



7 4 (e) 



1 + (50.0126265 . . .)- + (360.321 . . .)er 



(D.8) 
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By the PSLQ algorithm, this can be converted to the following expression 



/4(e) 



16 1 



1357T 2 e 2 
117301 



1 Z 434 n ( 1 ' 
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The first two numbers can be calculated analytically from the requirement of the multiplicative 
renormalizability. These agree with the numerical values computed as above. In the last number, 
the coefficients of n 2 , and ty 2 (\) can be fixed by the assumption that eg g is a fractional 
number. Then the remaining fraction agrees with the numerical value to six significant digits. 
The last integral /5(E) is found to be 



0.11851851851... 



1 + (21 .8926 ...)£+ (571 .233 . . .) y 



(D.10) 



This is equivalent to the following expression: 

16 



J 5 e 
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We again checked the first two numbers from the renormalization of the underlying field theory. 
The last number is again determined in an analogous manner to (ID .91) and agrees with the 
numerical value up to six significant digits. 
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